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Asymmetric  light  diffraction  by  pulsed  ultrasonic  waves 
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Low-MHz,  continuous  ultrasonic  waves  traveling  in  a  transparent  medium  cause  light  to  be 
diffracted  into  discrete  diffraction  orders  when  light  and  sound  propagation  directions  are  normal 
to  each  other.  When  pulsed  ultrasonic  waves  are  used  the  diffraction  orders  split  into  secondary 
orders  which  are  asymmetric  with  respect  to  the  central  diffraction  order.  This  splitting  is  derived 
and  a  general  expression  provided  for  the  intensity  as  a  function  of  the  ultrasonic  pulse  Fourier 
spectra.  Examples  are  provided  which  demonstrate  the  degree  of  asymmetry  for  an  exponential 
driving  pulse  and  the  convergence  to  the  classic  Raman-Nath  results  when  the  pulse  approaches  a 
continuous  wave. 

PACS  numbers:  43.35.Sx,  43.20.Bi 


INTRODUCTION 

Optical  probing  of  ultrasonic  waves  had  its  origins  in 
1932  through  the  independent  observations  by  Debye  and 
Sears'  in  Washington  and  Lucas  and  Biquard2  in  Paris  that 
an  ultrasonic  beam  in  a  liquid  acts  like  a  diffraction  grating 
when  illuminated  by  normally  incident  light.  In  1935  Ra¬ 
man  and  Nath1-5  explained  these  observations  by  treating 
the  ultrasonic  beam  as  a  moving  phase  grating.  Their  theory 
successfully  predicted  the  diffracted  light  angular  distribu¬ 
tion  and  the  relative  intensities  in  the  diffraction  orders  as  a 
function  of  sound  intensity,  optical  and  ultrasonic  wave¬ 
lengths,  and  ultrasonic  beam  thickness. 

Refinements  of  the  Raman-Nath  theory  have  been 
used  to  investigate  a  broad  spectrum  of  theoretical  and  ex¬ 
perimental  conditions  related  to  light  diffraction  by  ultra¬ 
sound.  This  has  ranged  from  the  theoretical  investigation  of 
light  diffraction  by  superposed  ultrasound  by  Murty6  to  the 
use  of  optical  probing  for  the  investigation  of  the  growth  of 
higher  harmonics  in  finite  amplitude  progressive  waves  by 
Zankel  and  Hiedemann7  and  Breazcale  and  Hiedemann.* 
Hargrove9  extended  the  Raman-Nath  theory  to  include  the 
prediction  of  the  diffraction  pattern  for  arbitrary  ultrasonic 
waves  as  illuminated  by  Gaussian  light  beams.  Zitter10  in 
turn  applied  this  extended  theory  to  light  diffraction  by  sym¬ 
metric  short  ultrasonic  pulses.  Other  extensions  have  includ¬ 
ed  the  examination  of  light  modulation  by  ultrasonic  waves 
in  the  presence  of  amplitude  optical  gratings  by  Calligari  et 
a/."-'3 


vides  a  direct  technique  for  the  relative  measurement  of  the 
ultrasonic  pulse  amplitude  spectra.  Second,  the  splitting  of 
the  intensities  in  the  first  Raman-Nath  diffraction  order  into 
satellite  orders  as  discussed  by  Zitter10  also  occurs  for  the 
central,  second,  third,  etc.,  diffraction  orders.  A  general 
expression  is  provided  for  the  calculation  of  these  intensities. 
Third,  the  intensity  distribution  within  the  diffraction  orders 
is  asymmetric,  i.e.,  the  intensity  in  the  positive  /nth  diffrac¬ 
tion  is  not  equal  to  the  intensity  in  the  negative  mth  order. 

I.  GENERAL  THEORY 

When  a  monochromatic  light  beam  of  frequency  /, 
wavelength  A,  and  width  21,  is  normally  incident  on  an  ultra¬ 
sonic  beam  of  diameter  D,  frequency/,,  and  wavelength 
the  amplitude  distribution  of  the  light  in  the  farfield  is  given 
by  the  diffraction  integral  as 

A  (9,t )  »  O'**  J"  (1) 

where 

C  =  the  normalization  constant 
a  =  Itrf 
k  =  2v/k 

6  =  the  farfield  angle 
( =  time 
w0  ■  2r/0 


Our  current  interest  is  prompted  by  the  recent  experi¬ 
mental  work  of  Hausler  et  a/.,14  which  demonstrated  the 
production  of  structured  diffraction  patterns  by  pulsed  ul¬ 
trasonic  waves  and  provided  qualitative  verification  of  the 
pulsed  diffraction  theory  presented  by  Zitter.10  This  paper 
expands  the  examination  of  Raman-Nath  diffraction  by 


*  Itr/Ag 
a  -  ~  *o* , 

with  v(a)  representing  the  phase  change  in  the  light  wave 
front  due  to  the  ultrasonic  wave.  For  a  continuous  wave  exci¬ 
tation 


pulsed  ultrasonic  waves  to  include  nonsymmetric  pulses  v(a)«  vsin(a),  (2) 

which  are  composed  of  modulated  sinusoidal  waves  of  wave- 
length  A,0  with  a  pulse  repetition  interval  A,. 

Thereare  three  primary  results.  First,  for  small  values  v«2 wpD/A.  (3) 

of  the  pulse  Raman-Nath  parameter,  optical  probing  pro-  defined  as  the  Raman-Nath  parameter  where  >i  is  the  maxi- 
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mum  variation  of  the  media  refractive  index.  For  an  arbi¬ 
trary  pulsed  ultrasonic  wave  with  functional  iomf{x,t,tp) 
where  tp  is  the  pulse  repetition  period,  v  can  be  expressed  in 
terms  of  the  Fourier  expansion  of  f(x,t,tp )  as 

«c 

Mo,)=  £  v»  J .  (4) 

n  -0 

with 

(op  =  2ir/tp 
ap  =  <opt  -  kpx 
v.  =anv 
kp  —  (op/v 

v  =  ultrasonic  velocity, 

where  an  and  <f>„  are  the  amplitude  and  phase  of  the  nth 
Fourier  component  of  the  pulse.  In  this  case  n  is  determined 
by  the  pulse  peak  amplitude. 

Substituting  Eq.  (4)  into  Eq.  (1)  and  moving  the  dc  com¬ 
ponent  of  v{ap)  outside  of  the  integral  yields: 

A  (0,t )  =  CefaVv°  j"' 

X  fl  e^^dx.  (5) 

#»  m  1 

Using  the  identity10 

i  (6) 

-  80 

the  product  expression  in  Eq.  (5)  can  be  rewritten  as 

fill-  2  V-™’’  (7) 

i»  —  I  m«  -  ® 

where  /,(e)  is  the  rth  order  Bessel  function  and 


r2*  ~  ® 

Xexp(  i{r,(m)di  +  “•  +  '«*„+•••)),  («) 

with  r,(m)  =  m  -  2r2 - nr„ - .  Substituting  Eq.  (7) 

into  Eq.  (5)  and  performing  the  integration  yields: 


with  A0  =  v0  sin  do  and 


(9) 


flm  =  (k  sin  0  -  mkp)l .  (10) 

Since  the  Vm  terms  are  complex,  Eq.  (9)  can  be  rewritten  as 


0-~  £  coa[(a»  +  m»,)r  +  40] 

-  V'"  sin[(a>  +  rruop)t  +  40] } 

+  /{ tf'jT  sin[(a>  +  mtop)t  +  40] 


+  #?«[(*  +  «•,)»  + 4.]  1)-^-=.,  (11) 

where  Re  and  Im  denote  the  real  and  imaginary  components 
of  tm.  The  light  intensity  is  then  given  by  the  square  of  the 
real  part  of  Eq.  (11). 

/(ftO-IRe^fti)}2.  (12) 

Toohtate  the  measurad  response  of  a  fret  photodetector 


it  is  necessary  to  take  the  time  average  intensity  over  an  in¬ 
terval  r  which  is  long  compared  to  the  period  of  the  incident 
light  and  at  the  same  time  is  short  compared  to  the  pulse 
repetition  rate,  i.e., 

7(0.*o)  =  —  P  +  l(0j)dt ,  a>—  >*Op  .  (13) 

T  Jh  T 

Normalizing  Eq.  (13)  to  the  peak  light  intensity  in  the  ab¬ 
sence  of  an  ultrasonic  wave,  i.e.,  0  =  0,  vm  =  0  for  m> 0, 
yields 


7  (0,to)  =  £  £ 

m  -  —  oo  n  *  — 


sin  fi„  sin  fl„ 

AT 


x({  +  d'Lmtf'im)cos(«  -  n)copt0 

-  I  *£■*?•  -  *'*'¥'im}sin(m  -  »K'o) . 

(14) 


which  is  independent  of  v0.  The  same  results  could  have  been 
obtained  by  incorporating  v0  into  the  unperturbed  media 
index  of  refraction.  Due  to  the  sin  x/x  term  in  Eq.  (14)  the 
light  intensity  will  exhibit  maxima  when  fl„  —  0,  i.e.,  at  the 
angles  0m  which  satisfy  the  relationship 


sin  6„=±m{AMp).  (15) 

Equation  ( 1 5)  is  the  same  as  the  result  obtained  in  continuous 
wave  diffraction  theory3  except  that  the  diffraction  angle  is 
determined  by  the  ratio  of  the  incident  light  wavelength  to 
the  pulse  spatial  duration.  Setting  n  =  m  —  1  and  evaluating 
the  product  (sin  fim  /(lm )  (sin  /II  „ )  at  the  angle  0  =  mA/ 
Ap  indicates  that  the  diffraction  orders  will  start  to  be  resol¬ 
vable  when  (2 nl/Ap)>  1.  Thus  as  the  pulse  repetition  rate 
increases  and  the  number  of  pulses  within  the  illumination 
interval  increases  (i.e.,  2irl/Ap>l),  (sin  )(sin  f2„/ 

/2„)~0  for  m#n,  and  the  diffraction  pattern  becomes  dis¬ 
crete,  i.e., 


l(0,t)  =  ?(0)=  £  ■  (16) 

where 


/.-(W+ow-iipj1.  ini 

with  Im  being  the  intensity  in  the  mth  diffraction  order. 

To  illustrate  the  behavior  of  the  diffracted  light  intensi¬ 
ty  as  the  ultrasonic  pulse  amplitude  increases,  *Pm  can  be 
factored  into  a  series  of  approximate  expressions  based  on 
successively  higher  Bessel  function  products,  i.e., 

A  \J \J  fJiJi,  etc. 


A.  First  ordsr  spproxlnutlon 

To  determine  how  the  diffraction  order  intensities  ini¬ 
tially  develop  we  can  examine  Eq.(  17)  in  the  limit  that  the  v„ 
are  small,  i.e.,  products  in  9m  of  order  S|(v,)/|(vk )  or  higher 
vanish.  Using  this  constraint  and  the  relationship 

m  »  r,  +  2/-,  +  -  +  nr.  +  - ,  (18) 

used  in  the  definition  of  fm,  we  find  that  for  m  *  0 

/.-In  (i9) 

I  f|na  I  I 

and  for  |m|>0. 
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i±m = I /±1(vm^±¥-  n  ^vjl1.  (20) 

I  I 

Since  /„  =  /_  m,  the  diffraction  pattern  is  symmetric  about 
the  central  order.  In  the  very  small  amplitude  limit  /0=sl.O 
and  I±m  ~vJ,/4.  Thus  the  intensities  in  the  individual  dif¬ 
fraction  orders  provide  a  relative  measure  of  the  ultrasonic 
pulse’s  Fourier  amplitude  spectra. 

.  _ _ I 


B.  Higher  order  epproxltwetlooe 

As  the  pulse  amplitude  increases,  the  JJx  and  J2  terms 
>n  *m  begin  to  contribute  and  the  intensity  in  the  mth  dif¬ 
fraction  order  ceases  to  be  solely  determined  by  y,(vm )  for 
m  >0.  For  the  central  diffraction  order  (m  =  0)  there  are  no 
additional  terms  due  to  the  constraint  imposed  by  Eq.  (18). 
However,  for  the  higher  diffraction  orders,  i.e.,  |m|  >0,  we 
find  that  the  intensity  is  given  by 


I±m=  n  Uvn)+J&mn\'±°*mn  II  Un) 

i»#m  k  yfcm/2 

+  2  k.lv^±,(v,_,)r±W’  +  *"-'  II 

« ■»  1  \  k  ,m  —  m  / 

n<m  —  *i 

+i( n  ^v*))|2, 

n  V  +  /  I 


(21) 


with  the  convention  that  all  terms  vanish  with  noninteger 
summation  indices  and  when  the  indices  of  the  v,  are  equal 
and  that  the  summations  and  II  products  are  valid  only  over 
positive  nonzero  index  values.  The  upper  sign  in  the  expres¬ 
sion  corresponds  to  the  positive  diffraction  orders  and  the 
lower  sign  corresponds  to  the  negative  orders. 

Since  J  _  „(x)  =  (  —  1  )V„  (jc),  the  second-order  terms  in 
(21),  J2(vm/2 )  and  the  two  J,J,  product  expressions,  exper¬ 
ience  a  1 80*  phase  shift  relative  to  Jx(vm)  in  going  from  +  m 
to  —  m.  If  these  terms  represent  a  net  increase  in  amplitude 
when  added  to  yv(vm ),  there  will  be  a  corresponding  net  de¬ 
crease  when  added  to  J  _  ,  lvm ).  This  causes  the  diffraction 
pattern  to  become  asymmetric. 

We  can  use  Eq.  (21)  to  illustrate  the  mechanisms  in¬ 
volved  in  the  growth  of  the  farfield  diffraction  pattern.  Based 
on  Eq.  (IS),  each  diffraction  order  corresponds  to  a  multiple 
of  the  pulse  repetition  rate  top.  The  first  term  in  Eq.  (21) 
represents  the  first  Raman-Nath  diffraction  order  for  the 
mth  Fourier  component  of  the  pulse.  The  second  term  repre- 
_ _ _ J 


I - - — “ 

sents  the  second  Raman-Nath  diffraction  order  of  the  m/ 
2th  pulse  Fourier  component.  The  J,J,  and  Jyf-i  summa¬ 
tions  correspond  to  the  mixing  of  the  contributions  of  the 
pulse  Fourier  components,  i.e.,  the  sum  and  difference  terms 
which  result  in  map.  The  resultant  farfield  diffraction  pat¬ 
tern  is  the  superposition  of  the  Raman-Nath  diffraction  pat¬ 
terns  of  the  pulse  Fourier  components  and  the  mixing  of 
these  components  as  constrained  by  Eq.  (18). 

In  the  continuous  wave  limit  all  Fourier  components  of 
v  vanish  except  for  the  v.  which  corresponds  to  the  pulse 
fundamental  frequency.  The  resulting  intensities 
/0  =  /oK);/±.  — i(v>i);  and  /±2*  =y|(vj  correspond 
to  the  zero,  first,  and  second  Raman-Nath  diffraction  or¬ 
ders. 

The  next  highest  approximation  adds  JtJPfirJrJ2>  *®d 
7,  terms  to  the  expression  for  the  intensity.  In  the  interest  of 
brevity  only  the  additions  to  Eq.  (19)  and  Eq.  (21)  are  shown 
with  the  phase  and  II  product  expressions  suppressed.  In 
each  case  the  additions  are  to  the  terms  inside  the  absolute 
value  signs.  The  central  order  is  given  by 


/„*  Eq.  (19)+  X  l*'-j(t',)y,(v2.)  +y2(v.)y_l(vto)) 

m  m  | 


2  •'-.Ky-.K-J 

\+y-,(vt/  x 

1  Vf-«  > 

1  )\ 

with  the  convention  that  all  terms  vanish  with  noninteger  summation  indices  and  when  the  indices  of  the  v,  are  equal  and  that 
the  summations  are  valid  only  over  positive  nonzero  index  values.  The  change  for  the  higher  orders  (|m|  >0)  is  given  by 


I±rn  -  Eq-  (21)+  X  /±j(V»)/±l(V-.-2.)  +  «/±3(V«/3)+  2  ;±j(V»^llV»— .1+  2  /T2(V»y±l(V2«  +  J 


+  2  (J±'W  2  /*i(v„yT,(r,-w_.)+/:fl(vf)  X  •'±i(v«y±i(vf+«.-.A 

f-l 1  A- I  A-l  I 

*  a<«  -m- A  ft<«  +  m  —  a  ' 


»<• 

«  —  I  m  ■  I 

y  +  xm-lf-a 


(23) 
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using  the  same  conventions  as  were  imposed  on  Eq.  .(21). 
Since  J _  „(x)  equals  ( —  l)Va(x),  the  /»  /j/„  and  /,/,/,  ex* 
pressions  maintain  a  constant  phase  relationship  with  Jx(vm) 
in  going  from  +  m  to  —  m.  As  a  result,  the  associated  in¬ 
crease  or  decrease  in  amplitude  for  the  +  m  and  —  m  orders 
is  the  same.  In  the  continuous  wave  limit  when  all  compo¬ 
nents  of  v  vanish  except  v.  the  diffracted  light  intensities 
become  /0  =  /oK);  /±„  =/|(v„);  I±U=J *nd 
I±Jm  =/}(»,).  These  correspond  to  the  zero,  first,  second, 
and  third  Raman-Nath  diffraction  orders.  Higher  order 
terms  will  impact  on  the  asymmetry  in  the  diffraction  pat¬ 
tern  only  if  the  sum  of  the  Bessel  function  orders  is  even. 


II.  ANALYTIC  RESULTS 

This  section  provides  examples  of  the  diffracted  light 
intensity  distribution  as  a  function  of  v  and  the  growth  and 
decay  of  asymmetries  within  the  diffraction  pattern.  For  il¬ 
lustration  the  exponentially  damped  sinusoid  given  by  Eq. 
(24)  is  used  as  the  pulse  time  history. 

t/(r)  =  ve~“*/2wk'  aiaadt  |A(f ) 

—  A  (/  —  2irAI/a»0))e,/4*' ,  (24) 

whereA  (f )  is  the  Heaviside  function,  (o0  =  2«/0and,  A,  and  k2 
determine  the  pulse  decay  and  repetition  rates.  In  this  form  a 
pulse  with  A,  =  10.0  and  k2  ~  30.0  has  an  e-folding  time  of 
10ro  and  a  repetition  rate  of  30r0  where  r0  =  l//0.  The  Four¬ 
ier  amplitudes  of  v  are  v„  =  vo„  withn.  =(<£  +  Aj),/2and 
e„  and  A.  given  by 


f  *~4M sin g.+e,  cose,} 

V  + 

Yn -e-^IJriny.  +  r,  cos  y,  ]  \  l/4*, 


(25) 


.  rin<„  -JcoseJ  +4 

‘•‘l - 7*7?. - 

g~4(r.  -^coay.)  +  A\ 

dl  +  ri  r  ' 


where  A  —  A2/A„  f,  =  2ir(A2  —  n),  and  y„  =  2v\k2  +  n). 
The  phase  »  given  by  arccos(c„/a„)  with  the  sign  of 
determined  by  the  quadrant  within  which  c,  and  b„  fall. 

Figure  1  illustrates  the  amplitude  and  phase  spectra  of 
Eq.  (24)  when  A,  =  10.0  and  k7  *  30.0.  For/,  =  3.0  MHz 
the  spectral  components  are  spaced  100  kHz  apart,  the  maxi¬ 
mum  amplitude  occurs  at  3.0  MHz,  and  the  components 
above  and  below  3.0  MHz  are  approximately  tr/2  out  of 
phase  with  a  slight  asymmetry  in  the  amplitude  spectra. 


A-  rlfMn  (MllUIOn  piiWnii 

Using  the  spectra  presented  in  Fig.  1  we  can  now  exa¬ 
mine  the  growth  of  the  diffraction  pattern  as  a  function  of  the 
Raman-Nath  parameter.  This  is  shown  in  Fig.  2  based  on 
Eqs.  (22)  snd  (23)  for  selected  values  of  v  with  an  arbitrary 
intensity  cutoff  set  at  0.1%,  i.e.,  IE-3,  on  the  abscissa.  In 
each  part  of  Fig.  2,  the  diffraction  orders  are  shown  with  the 
expected  relative  spatial  separation. 

Initially  when  v  *  0.5  [Fig.  2(a)],  the  central  order  has 


UMIin  Ml 

FIG.  I.  (a) Puhe amplitude ipectn for  v -i  1.0,  k,  =  10.0, and k2  «=  30.0. (b| 
Pulse  phase  spectra  for  v  «  1.0,  *,  -  10.0,  and  *,  —  30.0. 


no  satellite  orders  above  the  threshold  and  the  pattern  is 
almost  symmetric.  The  positive  and  negative  diffraction  or¬ 
ders  I  ±r>  through  /  ±  J,  locally  follow  the  asymmetry  which 
exists  in  the  amplitude  spectra  [Fig.  1  (a)].  As  the  pulse  ampli¬ 
tude  increases  the  second-order  sum  and  difference  terms 
previously  discussed  begin  to  introduce  asymmetries  into 
the  pattern.  At  v=  1.0  [Fig.  2(b)]  the  local  asymmetry  in  the 
satellite  orders  around  /_  *,  has  decreased  at  the  same  time 
that  the  asymmetry  in  the  satellite  orders  around  I+x  has 
started  to  increase.  As  the  second-order  terms  continue  to 
increase  in  magnitude  the  pattern  becomes  locally  asymme¬ 
tric  around  /_  *  [Fig.  2(c)],  the  asymmetry  between  the  posi¬ 
tive  and  negative  orders  becomes  more  pronounced,  and  the 
satellite  orders  begin  to  appear  locally  about  the  central  or¬ 
der.  The  asymmetry  in  7± ,  through  I ± ,  results  from  the 
J  _ terms  in  Eq.  (2 1 )  being  the  same  order  of  magnitude  as 
/,(vj.  As  the  pulse  amplitude  continues  to  increase  [Fig. 
2(d)j  the  local  asymmetry  around  the  central  order  decreases 
as  the  asymmetry  increases  around  I±  x.  At  v «  2.5  [Fig. 
2(e)],  the  number  of  satellite  orders  around  the  central  order 
has  increased  at  the  same  time  that  some  of  the  satellite  or¬ 
ders  about  1  ±yo  have  dropped  below  the  intensity  threshold. 
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The  contributions  of  the  approximations  in  Eqs.  (20h 
(23)  to  Use  ftrikkS  diffraction  pattern  presented  in  Fig.  2  are 
iBortrated  in  Fig.  3  a*  a  function  of  v.  The  tint  approxima¬ 
tion  Eq.  (20)  is  given  by  the  dotted  line,  the  second  approxi¬ 
mation  Eq.  (21 )  by  the  dashed  line,  and  the  tlurdapproxima- 
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F1Q.  2.  (a)  FarSeM  diffraction  pattern  for  v  *  0.5,  *,  =  10.0,  and  *2  =  50.0. 
(b)  FartMd  diffraction  pattern  for  v  «  1.0,  ft,  *  10.0.  and  k2  =  30.0.  (c)  Far- 
(letd  diffraction  pattern  for  v  —  I.S,  kt «  10.0,  and  t,  =  30.0.  (d)  Farfidd 
diffraction  pattern  for  v  =  2.0,  k ,  *  10.0,  and  *, »  30.0.  (e)  Farffdd  diffrac¬ 
tion  pattern  for  v  =  2.S,  k,  *  10.0,  and  k2  •  30.0. 


tion  Eq.  (23)  by  the  aolid  line.  For  /*,  Eqs.  (20)  and  (21)  are 
approximately  equal  and  Eq.  (23)  remits  in  a  small  correc¬ 
tion  to  the  intensity.  The  apparent  divergence  of  the  approxi¬ 
mations  for  I„  and  In  can  be  understood  by  ‘"■■nHirg  the 
relative  magnitudes  of  the  individual  contributions  to  the 
amplitude  of  the  diffracted  light  intensity.  For  v  an  2.0  and 

\H.Nsjghters.»  and  w.&  Meyer  uwrsaerto  light  ffthacWan  too 


m  as  29  the  leading  terms  in  the  series  approximations,  Eqs.  beyond  those  presented  in  Eq.  (23)  will  produce  second-or- 
(20)-(23),  excluding  the  phase  and  T1  product  expressions  are  der  corrections  to  the  intensity  over  the  range  of  interest.  For 
/,(vw);  J_  i(vt)/|(v30);  and  /_  i(v30)J,i(v2g)Jp1(vj1).  At  the  same  v>  1.0,  m  <  28  or  m  >  28  become  dominated  by  the  product 


time  the  validity  of  the  approximation  presented  in  Fig.  3  terms  that  correspond  to  the  frequency  associated  with  the 

can  be  gauged  by  examining  the  leading  terms  of  the  next  order.  For  v  above  3.0  higher  order  approximations 

higher  series  approximations,  /_  ,(v2)/_  ■(v30)/|(v29)/l(v3,)  JXJXJXJXJX;  etc.)  will  begin  to  produce  first-order 


and  /_,(v27)/1(v?g)/1(v29)/1(vj0)y_1(v31).  The  associated 
magnitudes  starting  with /((v^)  are:  1.4X10-1;  3.2xl0~3; 
3.3  X  10“3;  6.1  X  10~5;  and  2.5  X 10-5,  respectively. 

For  v  =  2.0  and  m  =  28,  the  leading  terms  in  the  series 
expansion  for  the  amplitude  are:  /,(v2g);  i(v2)/i(v30); 

^-i(v3oM,i(v27^i(vj,);  ,(v3y_  .(Vjol/^vwl/.tvj,);  and 

J-  i(v27^i(v2g)/i(v2,)/1(v30^_  ,(v32).  The  magnitudes,  start¬ 
ing  with  y,(v2g),  are:  7.8X10  2;  3.2x10  3;  2.3x10  3; 
6.1  x  10-5;  and  9.7 x  10  ‘,  respectively.  Based  on  the  ex¬ 
amination  of  /j,  and  /2g  it  is  apparent  that/,(vm )  for  decreas¬ 
ing  m  drops  in  magnitude  more  rapidly  than  the  leading 
terms  for  the  higher  approximations.  Also,  approximations 


181 


J.  Aooust  Sue.  Am.,  VW.  74,  No.  1 ,  Juty  1SW 


corrections  to  the  intensity. 

The  development  of  asymmetries  in  the  associated  neg¬ 
ative  orders  is  shown  in  Fig.  4.  As  expected  the  contributions 
of  the  /,/,  and  J2  terms  in  Eq  (21)  which  add  to  /29  (Fig.  3) 
result  in  an  equivalent  reduction  in  /_2,  (Fig.  4).  Concur¬ 
rently  the  higher  order  terms  in  Eq.  (23)  produce  a  reduction 
in  both  the  positive  and  negative  order  intensities. 


C.  Continuous  wav*  limit 

In  the  limit  that  the  pulse  approaches  a  continuous 
wave  the  satellite  orders  vanish  as  shown  in  Fig.  5  for  the 
central  order  (v  =  2.0).  As  kx  increases  the  pulse  amplitude 
spectra  approaches  vm  =  1.0  for  m  =  30  and  vm  =  0.0  for 
m  ^  30.  This  causes  / ,  ( vm )  to  decrease  more  rapidly  than  the 
;.,|v,)/,(v„  +  J  terms  locally  around  The  pattern  be¬ 
comes  symmetric  locally  about  the  central  order  and  then 
vanishes  as  the  central  order  intensity  approaches  J  J(2). 


PK}.4CbatritaUoaorBqs.(20H23)to/-Mthraueb/_na»atactioacfr 
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III.  EXPERIMENTAL  RESULTS 


To  provide  qualitative  verification  of  the  theory  pre¬ 
sented  in  Sec.  I,  an  Arenberg  pulser  was  used  to  excite  a  PZT 
transducer  with  a  short  pulse  (3.0-MHz  fundamental  fre¬ 
quency).  The  resulting  farfield  diffraction  pattern  from  a  se¬ 
quence  of  pulses  is  shown  in  Fig.  6  for  a  pulse  duration  of  ~4 
p s,  a  repetition  rate  of  —Him,  and  a  Raman-Nath  param¬ 
eter  of  — 3.0.  As  predicted  by  the  theory  the  splitting  occurs 
in  the  central  diffractioc  order  and  local  asymmetry  appears 
around /±M,  the  order  amocuted  with  the  pulse  fundamen¬ 
tal  frequency.  There  is  no  discernible  asymmetry  for  the 
satellite  orders  about  the  central  order  as  would  be  expected 
for  v~  3.0  due  to  the  dominance  of  the /(/_  |  terms  in  Eq. 
(21). 
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FIG.  6.  Typical  pulse  diffraction  pattern  for  3.0-MHz  pulses  with  4>is  dura¬ 
tion  and  32  jis  repetition  rate.  Top:  photograph  of  orders  up  ton  =  60,  with 
n  =  48  corresponding  to  location  where  the  firet  order  of  a  cw  diffraction 
pattern  would  appear.  Bottom:  enlargement  of  orders  around  n  =  48. 


IV.  SUMMARY 

An  analytic  expression  has  been  derived  for  the  Ra- 
man— Nath  diffraction  by  pulsed  ultrasonic  waves  which  can 
be  applied  to  any  general  pulse  shape  amenable  to  spectral 
analysis.  In  general  the  resulting  diffraction  pattern  is  asym¬ 
metric  around  the  central  order  with  the  degree  of  asymme¬ 
try  dependent  on  the  spectral  composition  of  the  pulse.  A 
limit**!  sot  of  examples  has  been  provided  which  demon¬ 
strates  the  degree  of  asymmetry  for  an  exponential  driving 
pulse  and  the  convergence  to  the  classic  Raman-Nath  re¬ 
sults  when  the  puke  approaches  a  continuous  wave. 

ACKNOWLEDGMENT 

The  research  presented  in  this  paper  has  been  per¬ 
formed  while  one  of  the  authors  (T.H.N.)  was  a  visiting 
member  of  the  Georgetown  University  Physics  Department. 


ip.  Debye  and  F.  W.  Sean,  Proc  Natl  Acad.  Sci.  USA  II 63  (1932). 

2R.  Lucas  snd  P.  Biqusrd,  “Noovdk»pcopd^optil|und»»ohde»etd** 
liquidc*  sotimis  i  faction  dc»  onde*  uhrasonores.”  C.  ft-  Acad.  Sci.  (Pari*) 
195, 121  (1932). 

5C.  V.Runan  and  N.S.N.  Nath,  “The  diffractioo  of  bght  by  high  frequen¬ 
cy  tound  wave*:  Part  I,“  Proc.  Ind-  Acad.  Set  B,  406-412  (1935). 

4C.  V.  Raman  and  N.  S.  N  Nath,  The  diffraction  of  light  by  sound  wave* 
of  high  frequency:  Part  II,"  Proc.  Ind.  Acad.  Sci.  H.  4)3-420(19?*'. 

5C  V.  Raman  and  N.  S.  N.  Nath, ‘The  diffraction  of  light  by  high  frequen¬ 
cy  foiiwl  waves:  Part  III,"  Proc.  Ind.  Acad.  Sci.  HI,  75—84  (1936). 

6J.  S.  Murty,  “Theoretical  investigation!  of  the  diffraction  of  light  by  super- 
posed  ultrasonic  waves,”  J.  Acoust  Soc.  Am.  26, 970-974  (1954). 

’K.  L.  Zaakri  and  E.  A.  Hiedemann,  “Diffraction  of  light  by  ultrasonic 
warn  progressing  with  finite  but  moderate  amplitudes  in  liquids,"  J. 
Acoust  Soc.  Am.  31, 44-54  (1959). 

*M.  A.  Breazeale  and  E.  A.  Hiedemann,  “Investigation  of  progressive  ul¬ 
trasonic  waves  by  tight  refraction,"  J.  Acoust.  Soc.  Am.  30,  751—756 
(1958). 

*L.  E.  Hargrove,  "Effects  of  ultrasonic  waves  on  Gaussian  light  beams  with 
comparable  to  ultraaonic  wavelength,”  I.  Acoust  Soc.  Am.  43, 
847-851  (1968). 

I0R.  N.  Zitter,  “Ultrasonic  diffraction  of  light  by  short  acoustic  pulses,"  J. 

Acoust  Soc.  Am.  43, 864-870  (1968). 

■  ip  Cslligaris,C  rvmipni  snd  PCiuti-  “Modulation  effects  of  progressive 
uhnsonic  a*»d«  on  s  light  beam  of  variable  width,”  Nuovo  Cunento  B  8, 
157-166(1972). 

,lF.  Calligaris,  P.  Gutt.sndl-  Gabridli,  “Extended  theory  of  Bfiht  modula¬ 
tion  in  thin -screen  diffraction  by  ultrasound,"  J.  Opt  Soc.  Am.  63,  287- 
292(1973). 

up.  CaUigsris,  P.  Ciuti,  and  I.  GsbrieUi,  “Light  modulation  by  ultrasound 
associated  with  an  amplitude  optical  grating.”  Acustics  27,  307—313 
(1972). 

UE  Hauler,  W.  G.  Mayer,  snd  M.  Schwartz,  “Light  diffraction  by  ultra¬ 
sonic  pulses,”  Acoust  Lett.  4, 180-184  (1981). 


fig  JlAmm'  9 MSga/.W' 


T.  H.  HOfipbort,  H  find  W.  ft.  iMyqr.  U»wr<ofi|l»  UK 


1 


Theoretical  prediction  of  a  hackscattering  maximum  at  Rayleigh 
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A  numerical  integration  method  for  the  description  of  acoustic  bounded  beams  is  used  to 
calculate  possible  backscattering  strength  from  a  smooth  liquid-solid  interface.  It  is  shown  that 
the  backscattering  strength  is  masimum  for  Rayleigh  angle  incidence.  The  influence  of  beam 
shape  and  beam  width  on  the  backscattering  strength  near  the  maximum  is  demonstrated. 

PACS  numbers:  43.33.Pt,  68.23.  +j,  43.20.Bi,  43.20.Fn 


INTRODUCTION 

Recently,  an  interesting  phenomenon  associated  with 
the  backscattering  of  a  bounded  acoustic  beam  from  a 
smooth  liquid-solid  interface  has  been  reported. 1  This  phen¬ 
omenon  involves  the  existence  of  a  relative  maximum  at  the 
Rayleigh  angle  incidence  when  the  backscattering  strength 
is  measured  for  various  incident  angles.  First,  de  Billy,  Ad¬ 
ler,  and  Quentin1  reported  on  this  maximum  of  backscatter¬ 
ing  strength  for  a  water-stainless  steel  and  a  water-copper 
interface.  The  same  authors  later  considered  a  more  general 
scattering  geometry  and  observed2  a  relative  maximum  in 
the  scattering  strength  at  both  backward-scattered  and  for¬ 
ward-scattered  Rayleigh  angles  for  a  general  incident  angle. 
Indeed,  these  observations,  described  as  conical  reflection, 
have  been  qualitatively  reported3  earlier  by  Diachok  and 
Mayer  for  the  specific  case  of  Rayleigh  angle  incidence. 

These  experimental  results  have  prompted  several  theo¬ 
retical  investigations  attempting  to  explain  these  phenom¬ 
ena.  Within  the  context  of  spectral  representation  of  acoustic 
bounded  beams,  Norris4  extended  the  Bertoni  and  Tamir 
analytical  procedure,5  developed  for  forward  reflection,  to 
the  backscattering  case.  In  light  of  the  fact  that  surface 
roughness  should  play  a  significant  role  in  backscattering 
measurements,  de  Billy  and  Quentin6  employed  the  poten¬ 
tial  method  as  described  by  Welton7  to  evaluate  the  acoustic 
field  of  a  bounded  beam  scattered  from  a  rough  liquid-solid 
interface.  Comparing  their  theoretical  results  with  back- 
scattering  measurements  on  samples  with  rough  interfaces, 
de  Billy  and  Quentin  showed  good  agreement  for  incident 
angles  of  4(T  or  less  except  that  their  theory  did  not  account 
for  the  relative  maximum  observed  near  the  Rayleigh  critical 
angle. 

The  above  theoretical  efforts4,6  appear  to  confirm  the 
belief  that  the  relative  maximum  of  backscattering  strength 
at  Rayleigh  angle  incidence  must  be  associated  with  the 
boundedness  of  the  acoustic  source  as  well  at  the  resonant 
surface  propagation  mode  corresponding  to  the  Rayleigh 
critical  angle.  Under  this  motivation  the  present  study  in¬ 
tends  to  investigate  this  backscattering  phenomenon  of  a 
smooth  liquid  solid  Interface,  using  an  exact  numerical  inte- 

scribed  by  the  spectral  representation.  This  numerical  algo¬ 


rithm  was  initially  developed  by  Ngoc  and  Mayer"  to  de¬ 
scribe  acoustic  bounded  beams  reflected  from  or  transmitted 
through  layered  media.  As  will  be  shown  in  the  following 
calculation  results,  this  theoretical  analysis  is  able  to  de¬ 
scribe  the  general  features  of  backscattering  strength  as  a 
function  of  incident  angle,  including  the  relative  maximum 
at  the  Rayleigh  critical  angle. 

I.  THEORY  AND  COMPUTATIONAL  RESULTS 

Following  the  formulation  developed  by  Ngoc  and 
Mayer,"  the  spectral  representation  is  again  invoked  to  de¬ 
scribe  the  scattered  field  of  an  acoustic  bounded  beam  of 
frequency  to  inciden  t  ">nto  a  smooth  liquid-solid  interface  at 
an  incident  angle  0,, 

U\x*)  =  {lirr'^R(k,)V(kx) 

Xexp[i(x*x  +**,)]  dkx.  (1) 

The  symbols  used  in  above  equation  are  defined  as  follows: 

U  (*,*)  is  the  acoustic  scattered  field  in  the  (x,z)  plane  assum¬ 
ing  uniformity  in  they  dimension;  k  is  the  wave  vector  in  the 
liquid;  kx  and  kt  are  the  x  and  z  components  of  the  k  vector 
related  by  kt  =(**  —  kj),/2.  The  quantity  R  [kx )  denotes 
the  plane-wave  reflection  coefficient  at  a  liquid-solid  inter¬ 
face  (e.g.,  see  Ref.  8)  and  K(kJ  is  the  Fourier  transform  of  the 
incident  acoustic  beam  given  at  the  interface,  describing  the  I 

complex  amplitude  of  the  constituent  plane  waves  that  form  | 

the  incident  beam.  ( 

In  calculating  the  scattered  sound  field  the  integral  of 
Eq.  (1)  is  interpreted  as  describing  the  scattered  field  repre-  j 
sen  ted  by  a  superposition  of  plane  waves  having  the  complex 
amplitude  R  ).  Integrating  from  —k  tok  implies  | 

that  contributions  to  die  scattered  field  can  come  from  plane 
waves  reflected  in  all  directions  from  this  interface.  To  be 
compatible  with  most  experimental  arrangements,  trans-  i 
ducers  of  finite  size  for  both  acoustic  source  and  detector  will  I 

be  taken  into  account  in  the  present  study.  In  addition,  one  J 
would  also  like  to  be  able  to  account  for  the  specific  response 
pattern  of  the  detector  in  use.  The  mattering  geometry 
adopted  in  the  following  computations  is  Jhwtratad  in  Fig. 

1- The  source  and  detecting  transducers,  both  havingawidth 
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FIG  1.  Schematic  dacripdoc  of  the  nattering  geometry  under  comidert- 


of  2 w,  are  rotated  on  semicircular  paths  of  radius  d.  The 
incident  and  detecting  angles  are  defined  as  0,  and  0d,  re¬ 
spectively.  The  positive  direction  of  these  angles  is  indicated 
by  the  underlying  arrows  and  therefore  detection  of  the 
backscattered  sound  would  assume  a  negative  value  for  9d . 
For  the  backscattering  geometry,  0d  —  —0,. 

In  the  backscattering  calculations  to  be  presented,  the 
incident  beam  profile  is  taken  to  be  a  finite  Gaussian  one  and 
the  detecting  transducer  is  assumed  to  have  a  flat  response. 
The  calculations  are  carried  out  according  to  the  following 
procedure: 

(a)  For  a  given  0„  the  incident  beam  profile  is  first  pro¬ 
jected  onto  the  r  =  0  interface; 

(b)  The  projected  incident  profile  is  then  Fourier  trans¬ 
formed  to  yield  V(kx ); 

(c)  Interaction  of  individual  plane  waves  of  complex 
amplitude  V(kx )  with  the  interface  is  described  by  the  pro¬ 
duct  R  (kx )  V(kx ),  which  is  now  entered  into  the  integrand  of 

Eq-(l); 

(d)  The  scattered  acoustic  field  at  spatial  position  (x,z) 
can  now  be  evaluated  by  an  exact  numerical  integration  of 
Eq.  (1)  across  the  surface  of  the  detecting  transducer,  which 
can  be  represented  by 

z  =  tan  9d\d  —  x) 
for 

d  tin  6d  —  to  cos  9d  <jc<d  sin  ed+u>  cos  e4\  (2) 

(e)  The  backscattering  strength  for  a  given  9d  =  —  0,  is 
finally  determined  from  the  scattered  field  calculated  across 
the  detecting  surface  in  the  preceding  step  by  averaging  them 
according  to  the  selected  response  of  the  detector. 

Figure  2  presents  the  first  set  of  computational  results 
showing  backscattering  strength  as  a  function  of  incident 
angle  far  three  values  ofbesmwidth,  2w  *  3, 10,  and  20  mm, 
with  an  acoustic  source  of  2  MHz.  These  computations  are 
done  for  a  water-stainless  steel  interface  with  the  radial  dis¬ 
tance  d  being  30  mm.  The  backscattering  strength  curves 
shown  in  Fig.  2  exhibit  a  steady  decrease  as  the  incident 
angte  becomes  larger,  with  adistinct  relative  maximum  posi¬ 
tioned  approximately  at  the  Rayleigh  critical  angle  whidi  is 
calculated  to  be  30.65*  in  this  case.  Near  the  relative  maxi¬ 
mal  one  observes  that  the  maximum  peak  increases  with 
bsaanwidth  and  the  slope  of  the  maximum  is  broadened  at 


10.0  20.0  30.0  40.0 

INCIDENT  ANGLE 


FIG.  2.  Variation  of  baducattcring  strength  with  beamwidth  for  a  Suite 
Gaussian  beam  of  2  MHz  incident  on  a  water-stainless  steel  interface 

Next,  the  bcamwidth  is  kept  constant  at  the  value  of  10 
mm  but  the  shape  of  the  incident  Gaussian  profile  is  varied. 
This  is  achieved  by  introducing  a  parameter  a  in  the 
expression  describing  a  Gaussian  profile, 

tftac  W  =  exp  -  [{x/aw)2].  (3) 

Setting  a  =  1  corresponds  to  the  incident  profile  used  in 
the  previous  set  of  calculations.  In  the  present  set  of  calcula¬ 
tions,  the  incident  beam  generated  by  setting  a  =  0.5  has  a 
narrower  profile  while  a  —  2.0  results  in  a  broader  profile. 
These  incident  profiles  are  illustrated  in  Fig.  3. 

Again,  the  above  computational  procedure  is  applied  to 
calculate  the  backscattering  strength  versus  the  incident  an¬ 
gle  for  the  three  incident  profiles  shown  in  Fig.  3.  The  results 
are  presented  in  Fig.  4  where  the  steady  decrease,  inter¬ 
rupted  by  a  relative  maximum  near  the  Rayleigh  angle,  is 
observed  for  all  of  these  profiles.  The  most  noticeable  differ¬ 
ence  among  them  is  that  the  relative  maximum  is  consider- 
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FIG.  4.  Variation  of  backscattering  strength  due  to  chancel  in  shape  of  the 
incident  profile.  The  experimental  curve  ii  redrawn  from  Ref.  1  for  the  case 
of  a  Gaussian  beam  of  4.S  MHz  and  12-mm  beam  width,  incident  on  a  ws- 
ter -stainless  steel  interface  with  other  parameters  not  specified. 


ably  more  pronounced  when  the  incident  beam  profile  be¬ 
comes  narrower.  A  typical  experimentally  obtained  curve  of 
backscattering  strength  is  redrawn  from  Ref.  1  and  is  includ¬ 
ed  in  Fig.  4  for  comparison. 

II.  CONCLUSION 

The  above  computations  establish  that  a  theoretical  de¬ 
scription  based  on  spectral  analysis  successfully  predicts  a 
relative  backscattering  maximum  at  the  Rayleigh  critical 
angle  and  other  general  features  of  the  backscattered  sound 
field  observed  by  recent  experimental  measurements.  A 


preliminary  application  of  this  theoretical  approach  has 
identified  two  key  parameters  that  would  affect  the  back- 
scattering  strength  as  a  function  of  the  incident  angle. 
These  two  parameters  are  beam  width  and  profile  shape  of 
the  incident  beam.  It  is  expected  from  this  formulation  that 
other  parameters  including  response  pattern  of  the  detec¬ 
tor  and  distance  of  source  and/or  detector  from  the  inter¬ 
face  would  also  strongly  influence  the  backscattering 
strength. 
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Abstract  -The  reflection  of  non-Gaewiaa  ultrasonic  beams,  expressed 
in  terms  of  nth  ordered  polynomials,  is  discerned  and  it  is  shown  that  at 
critical  kmgitiidiiial  and  RayM0i  aa0es  of  incidence  deviations  from 
the  Gaussian  profile  case  occur,  varying  with  the  polynomial  representa¬ 
tion  chosen.  Examples  are  given  for  water-Plexiglas  and  water-stainless 
steel  flat  interfaces. 

Introduction 

Nonspecular  reflection  and  transmission  effects  of  a  finite 
sound  beam  have  been  investigated  1 1 1— 1 7 J  for  various  types 
of  layered  media.  Bertoni  and  Tamir's  analysis  1 1 )  described 
these  effects  at  the  Rayleigh  critical  angle  for  a  liquid-solid 
interface.  Pitts  era/.  (2)  extended  Bertoni  and  Tamir's  ap¬ 
proach  to  another  structure  in  which  a  solid  plate  is  immersed 
in  a  liquid  medium.  This  approach  again  was  able  to  describe 
nonspecular  effects  at  the  critical  plate-mode  angles  associated 
with  Lamb  waves  of  the  solid  plate.  The  problem  was  general¬ 
ized  to  incorporate  sound  attenuation  in  the  media  and  was 
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solved  numerically  by  Ngoc  and  Mayer  |3]-[6) .  This  led  to  a 
description  of  nonspecular  reflection  and  transmission  effects 
at  all  incident  angles  where  these  effects  may  occur.  The 
numerical  approach  enabled  the  discovery  of  nonspecular 
effects  at  the  longitudinal  critical  angle  (3)  of  the  liquid-solid 
structure  and  at  the  between-mode  angles  (5)  of  the  liquid- 
solid-liquid  structure  in  addition  to  those  taking  place  at  the 
Rayleigh  angle  and  the  critical  plate-mode  angles  of  the  respec¬ 
tive  structures.  It  was  also  able  to  describe  the  variation  of 
nonspecular  features  as  the  incident  beam  was  steered  away 
from  the  critical  angles. 

The  results  of  these  studies  provided  an  adequate  under¬ 
standing  of  the  physical  processes  that  underlie  the  nonspec¬ 
ular  effects  taking  place  upon  reflection  from  or  transmission 
through  s  flat  interface.  However,  a  major  deficiency  of  the 
above  theoretical  models  that  limits  the  application  of  these 
effects  to  practical  problems  rests  with  the  fact  that  these 
analyses  are  all  based  on  the  assumption  of  an  incident  beam 
having  a  Gaussian  intensity  distribution. 

In  most  practical  situations  where  one  finds  a  sound  beam 
incident  at  an  interface,  the  incident  beam  profile  would  in¬ 
variably  be  non-Gaussian.  This  is  due  to  the  fact  that  all  acous¬ 
tic  beams  traveling  in  any  kind  of  medium  must  be  subjected 
to  several  mechanisms  that  will  distort  the  beam  patterns  (8) , 
(9) .  Such  mechanisms  include,  for  example,  nonlineer  inter¬ 
action,  absorption,  geometrical  diffraction,  and  dispersion  of 
the  medium  itself.  These  mechanisms  would  sometimes  change 
the  beam  pattern  produced  by  the  transducer  by  enhancing 
or  reducing  the  side  lobes.  It  is  therefore  quite  desirable  to  be 
able  to  understand  how  an  incident  beam  with  side  lobes  be- 
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haves  after  interacting  with  an  interface.  In  addition,  one 
would  also  like  to  know  if  sidelobes  are  generated  as  part  of 
nonspecular  reflection  or  transmission  of  a  non-Gaussian  inci¬ 
dent  beam.  As  a  result  one  can  distinguish  the  changes  in 
beam  profile  due  to  the  nonspecular  effects  from  those  pro¬ 
duced  by  the  above-mentioned  mechanisms  associated  with 
the  medium. 

In  general,  the  current  theoretical  models  describing  non¬ 
specular  effects  would  become  more  useful  if  they  allowed  for 
treatment  of  incident  beams  of  an  arbitrary  intensity  distribu¬ 
tion.  In  this  paper,  the  numerical  approach  developed  for 
Gaussian  incident  beam  will  be  modified  to  investigate  the 
nonspecular  effects  for  incident  beams  having  an  nth-order 
polynomial  intensity  distribution.  The  choice  of  the  poly¬ 
nomial  distribution  was  motivated  by  the  consideration  that 
polynomial  distributions  can  approximate  very  well  many 
transducers  in  current  use;  furthermore,  they  can  be  adapted 
to  be  used  as  basic  components  to  construct  an  arbitrary  in¬ 
tensity  distribution.  In  the  following,  modifications  to  the 
theoretical  framework  will  be  illustrated  and  sample  nonspec¬ 
ular  reflected  beam  profiles  computed  from  the  modified 
theory  will  be  presented. 

Theoretical  Formulation 

Consider  a  finite  beam  bounded  in  the  (x,  2)  plane  and  uni¬ 
form  in  the  >> -dimension,  of  angular  frequency  co  =  2itf  and 
beam  width  2w.  A  spectral  representation  of  such  a  bounded 
beam  treats  it  as  a  superposition  of  an  infinite  number  of  plane 
waves  having  the  same  frequency  but  different  amplitude  and 
incident  at  slightly  different  angles  about  a  central  direction 
indicated  by  0t.  Extending  the  principle  of  spectral  represen¬ 
tation  to  a  reflected  or  transmitted  bounded  beam  one  can 
describe  the  reflected  or  transmitted  field  distribution  by 
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INCIDENT  PROFILES 


Fig.  1 .  Profile  of  incident  beam  represented  by  zero-ordered  poly¬ 
nomial  [step  function)  (long  dashes);  2nd-ordered  polynomial  (short 
dashes);  4th-ordered  polynomial  (dash-dot);  6th-ordered  polynomial 
(dash-double  dot);  and  a  Gaussian  distribution  (solid  line). 

where  the  function  erf(x)  is  defined  by 

erf(jr)  =  2/(ff)1'1  f  exp  [-I1)  dr.  (6) 

do 

Combination  of  (3)  and  (5)  gives  the  normalization  constants 

AJn  up  to  the  sixth  order  as 

A0  *0.7734, 

A 2  *  1.0590, 

At  =  1.2132, 


•  exp  \i(xkx  +Z**))  dkx. 


(1)  A6  =  1.3243. 


(7) 


where  f\kx)  is  the  plane  wave  reflection  or  transmission  co¬ 
efficient  for  the  layered  structure  under  consideration.  As 
before  [5] ,  V(kx)  is  the  Fourier  transform  of  the  incident 
field  and  w0  =*  w/cos  8t  and  kt~  k  sin  0j. 

The  present  study  investigates  an  incident  beam  described 
by  an  nth-order  polynomial  and  compares  the  computational 
results  with  those  for  a  Gaussian  incident  beam.  Since  most 
sound  beams  in  practice  are  symmetric  only  the  even  orders  of 
polynomials  are  considered.  The  specific  form  of  the  poly¬ 
nomials  used  to  model  the  incident  beam  in  this  study  are 
taken  to  be 


Mjsll-Wwo)1]", 

|0, 


-w0  <  x  <  w0; 
elsewhere, 


(2) 


where  AiH  is  the  normalization  constant.  In  order  that  re¬ 
flected  beam  profiles  calculated  for  a  polynomial-type  incident 
beam  can  be  compared  to  those  for  Gaussian  beams,  A  2„  will 
be  normalized  according  to 


A\h  f"°  II  *  ix/w0)i]1"dx  -  f  °  I  Iff  I*  dx,  (3) 

*-»o  •'-Wo 

where  Uf(x,  0)  is  the  Gaussian  incident  distribution  given  by 
Vt{x,  0)  * exp  I-  (x/w0)1  +  fcr*|| .  (4) 

Evaluation  of  the  right-hand  side  of  (3)  yields 


\Ut\*  dx  ■  (ir/2)1^*  w0  erf(>/2), 


(5) 


The  normalized  intensity  profiles  of  the  even-ordered 
polynomial  incident  beams  and  the  Gaussian  intensity  profile 
are  presented  in  Fig.  1  for  comparison.  It  is  noted  that  the 
analysis  of  the  general  even-ordered  polynomial  beams  also 
covers  the  special  case  of  a  piston  source,  where  n  *  0. 

In  order  to  evaluate  the  reflected  or  transmitted  sound  field 
from  ( 1 ),  one  needs  to  determine  V(kx )  from  the  polynomial 
incident  profile  given  by  (2).  This  results  in 
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Fig.  2.  Rayleigh  angle  reflection  from  a  water-stainless  steel  flat  inter¬ 
face  for  a  Gaussian  incident  beam  (long  dash-dot);  zero-ordered  poly¬ 
nomial  (short  dashes);  2nd-ordered  polynomial  (short  dashes-dots); 
4th-ordered  polynomial  (short  dasheadouble  dots);  6th-ordered 
polynomial  (long  dashes). 

Evaluation  of  the  integral  in  ( I )  to  determine  the  reflected 
or  transmitted  sound  field  is  now  possible  since  Vlkx)  has 
been  determined  and  the  plane  wave  reflection  or  transmission 
coefficient  P{kx)  for  several  layered  media  is  well  established 

131. 

Computational  Results 

The  nonspecular  effects  for  an  incident  beam  having  an  even- 
ordered  polynomial  profile  are  illustrated  here  for  the  simple 
case  of  a  bounded  beam  reflected  from  a  liquid-solid  interface. 
Calculations  are  performed  for  the  Rayleigh  critical  angle  of  a 
water-stainless  steel  interface  and  for  the  longitudinal  critical 
angle  of  a  water-Plexiglas  structure. 

In  Fig.  2,  the  nonspecularty  reflected  beam  profile  at  the 
Rayleigh  angle  is  presented  for  both  types  of  incident  profiles. 
In  the  well-known  case  of  a  Gaussian  profile,  the  reflected  pro¬ 
file  shows  the  familiar  two  peaks  with  a  trailing  field  on  both 
sides  of  the  beam.  However,  reflected  profiles  associated  with 
the  even-ordered  polynomial  incident  beam  all  exhibit  distinct 
sidelobes  in  place  of  the  trailing  field.  Among  profiles  of  dif¬ 
ferent  orders,  the  profiles  of  lower  orders  have  a  higher  reflec¬ 
tion  intensity  level  as  measured  by  the  magnitude  of  the  two 
principal  peaks.  In  general,  one  observes  that  the  polynomial 
description  of  the  incident  beam  does  not  change  the  funda¬ 
mental  nonspecular  features  other  than  the  appearance  of  the 
sidelobes. 

For  the  water-Plexiglas  interface,  nonspecular  effects  are 
investigated  for  the  longitudinal  critical  angle.  Results  of  an 
earlier  study  [3J  for  a  Gaussian  incident  beam  only  showed 
the  nonspecular  feature  of  a  lateral  displacement.  For  the  case 
of  polynomial  incident  beams  the  nonspecularly  reflected 
beams  for  all  orders  considered  are  again  characterized  by  a 
lateral  displacement  of  approximately  the  same  magnitude. 
Sidelobes,  although  discernable  but  very  small,  are  also  present 
on  both  sides  of  the  main  peak.  These  results  are  illustrated  in 
Fig.  3.  The  appearance  of  sidelobes  in  both  sets  of  calcula¬ 
tions  is  definitely  associated  with  the  choice  of  an  nth-order ed 
polynomial  for  the  incident  beam  distribution.  For  such  a 
distribution,  decomposition  of  the  incident  beam  into  plane 


Fig.  3.  Critical  longitudinal  angle  reflection  from  a  water-Plexiglas  flat 
interface  for  a  Gaussian  incident  beam  (long  dashes-dots);  zero- 
ordered  polynomial  (short  dashes);  2nd-ordered  polynomial  (short 
dashes-dots);  4th-ordered  polynomial  (long  dashes). 

waves  as  represented  by  its  Fourier  transform  Fjn(a)  shows 
contributions  from  plane  waves  deviating  considerably  from 
the  central  direction  are  quite  significant  and  vary  according 
to  sin  a/a  or  cos  a/a.  As  a  result,  the  trailing  field  which 
always  exists  in  the  case  of  a  Gaussian  incident  profile  is  no 
longer  present  and  sidelobes  are  introduced  into  the  reflected 
profiles. 

Conclusion 

The  description  of  a  finite  incident  beam  by  an  even-ordered 
polynomial  has  been  shown  to  lead  to  a  new  nonspecular  fea¬ 
ture  in  the  reflected  beam  profile.  The  computational  results 
predict  that  upon  reflection  at  the  critical  angles  associated 
with  layered  media,  sidelobes  will  appear  in  the  reflected 
profile  if  the  incident  beam  has  a  polynomial  distribution 
instead  of  a  Gaussian  one.  The  present  study  can  form  a  basis 
for  investigations  of  nonspecular  reflection  or  transmission  of 
an  incident  beam  having  an  arbitrary  intensity  distribution 
which  can  be  modeled  through  polynomial  fitting. 
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